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Chapitre 5 

Intégration et dérivation numérique 



5.1 Introduction 

5.1.1 Objectif 

Le but est l'évaluation numérique de l'intégrale définie d'une fonction. Lorsque l'intégrale de la 
fonction / n'est pas calculable ( par exemple f Q e~ x dx) on verra que l'on peut l'approximer de 
diverses façons. Des fois même, connaissant la primitive F de / il est plus pratique d'utiliser ces 
méthodes d'approximation que de procéder à une évalution directe en fonction de F. 
Dans ce chapitre, nous allons présenter la principe général du calcul approché d'intégrale que nous 
appliquons dans le cas de la méthode des rectangles, trapèzes et Simpson. Nous présentons aussi 
la théorie des formules de quadrature. La dernière partie sera consacrée a la dérivation numérique. 

5.1.2 Rappels: 

i. Si / est continue sur [a, b], alors / admet une primitive donc fj: f(t)dt existe pour tous 
a,0e [a,6], 

ii. Soient / : }a,b[ — > IR et F: [a, b] — > IR : F est une primitive de / sur [a, b] si F est continue 
sur [a, b], si F est dérivable sur ]a,b\ alors F'(x) = f{%), quelque soit x ç]a.b[. 

iii. Si F est une primitive de / alors G = F + cte est aussi une primitive de /. 

vi. a-, f3 G [a, b] , si F et G sont deux primitives de / on a: 

/ = F(j3) - F(a) = G(J3) - G(q-) 

/ est indépendant de F et ne dépend que de a, ft et /. 

vii. - Oi.fi G [a, b], si F et G sont deux primitives de / on a: 

/ = F(f3) - F(a) = G(J3) - G(a) = f f(t)dt 

J a 

I est indépendant de F et ne dépend que de a, j3 et / 

viii. Si t G [a, 6], la fonction F de [a, b] dans IR définie par F(t) = f f(x)dx est une primitive de 
/: telle que F(a) = et F'(t) = f(t) pour tout t e [a, 6]. 

5.1.3 Interprétation géométrique de l'intégrale 

Soit une fonction numérique /: [a, b] — > IR continue, 

a. si / > sur [a, b], f a f(t)dt est la surface hachurée de la figure qui suit: 

b. Si / change de signe sur [a, b], dans ce cas on considérera la somme algébrique. 
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4- Intégration et dérivation numérique 



y=f(x) 



Figure 5.1: / de signe constant 



y=f(x) 



Figure 5.2: / change de signe. 



5.1.4 Principe du calcul approché: 

l'interprétation géométrique de l'intégrale suggère de décomposer la surface sous le graphe de la 
fonction en éléments d'aires simples (rectangles, triangles, trapèzes ...) moyennant certaines ap- 
proximations au voisinage de la courbe. 

L'idée est donc de subdiviser le segment [a, b] en plusieurs sous-segments [a^ aj+i] (i = 0, . . . , ri — 
l,ao = a,a n = b). 



/ f(t)dt = J2 / /(*)* 



on est donc ramené à évaluer ff i+1 f(t)dt sur un petit segnrent [a*, Uî+i] où / ne varie pas trop. 

Ce qui nous amène à proposer la formule approximative suivante sur un petit intervalle [ce, fi] 
de [a,b]: 

(5.1) 



/' f(t)dt = I aj3 if) + R a ^(f) 

Ja 



lafîif) est l'approximation proposée, R a ^(f) est un reste: "erreur de méthode" (5.1) est dite 
formule simple ou élémentaire, la formule (5.1) appliquée sur chaque segment [ctj, ai+i] nous donne 
la formule approximative Composée: 



/ f(t)dt = J2 4„a 1+1 (/) + E «*,*+>(/) 



Dans tout ce qui suit nous considérons la subdivision suivante: 

ai = a -\- ih\ 2 = 0,..., n, clq = a, ai = a + ih. a n = b avec h ■ 



b — a 
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y=f(x) 




Figure 5.3: Subdivision de [a, b] et approximation de f a f(x)dx par la somme des surfaces hachurées. 
cette subdivision est appelée la subdivision régulière de pas h 

5.2 Formule des rectangles ( avec point milieu ): 

5.2.1 Formule élémentaire: 

On remplace le graphe de la fonction / entre [a,/3] par une droite y = f{"f), où 7 = ■—^ est le 
point milieu de [a,/3]. (voir figure ci dessous) 





y=f(x) 






f(p) 






-"""Cf 


f(Y) 








f(cc) 






( 










a 


Y 


P 



Figure 5.4: Mj'ethode des réctagles avec point milieu. 



• 'a 



Estimation du reste R a ^{f): 

Supposons que / a une dérivée seconde continue ( F aura donc une dérivée tierce continue ). La 

fomule de Taylor à Tordre 3 entre 7 et (3 appliquée à F: 



HP) -F{l) = {l3-l)F'{l) + \{p-lfF"( 1 ) + ^ [ {li-afF'"{Q 
= |/(7) + y/'(7) + g/"(0; €€ 1)9,7] 



(5.2) 
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Formule de Taylor entre a, 7: 

F(a)-F(y) = -|/( 7 ) + Ç/'( 7 ) - ^/"W V € [0,7] (5-3) 

En retranchant (5.2) de (5.3) on obtient la relation: 

j" est une fonction continue, d'après le théorème des valeurs intermédiaires il existe 6 ée]ï7,Ç[ telle 
que 

\{f"{n) + /"(«)) = /"W 

La formule d' élémentaire devient donc: 

[\f(t)dt = hf( 7 ) + '^f"(0); a<0<0 (5.4) 

5.2.2 Formule composée: 

Supposons que: /: [a.b] — » ïïl a une dérivée seconde /" continue, considérons la subdivision 
(ai)i-o,i....7î de pas /t = ^p, sur chaque intervalle [a^Qi-i-i] on a la formule élémentaire: 



/ /(£)<& = hf("fi) + — / "(0i) ; ^ <9 { < a l+1 avec 7, ■■ 



a i + a i+l 



d'où la formule composée suivante 

rb 



fb "- 1 ]yi "- 1 1 

/ f(t)dt = hJ2 f(l,) + n-J2 -f"( 8 ') 
Or, /" est continue, d'après le théorème des valeurs intermédiaires, il existe 9 €E [a, b] tel que: 

" i=0 



en tenant compte du fait que nh = (b — a), on obtient comme formule composée: 

«=n-l h 2 

■A ' 9 '" ) " 
i=0 



/*/(t)* = ^^/(^ii) +g(6-a)/"(0) 



£55 B " (/) 

graphiquement on a le schéma suivant: 

5.2.3 Propriétés. 

a. Si /" = sur [a, b] c-à-d / est un polynôme de degré un ou moins, il n'y a pas d'erreur de 

méthode. 

b. Si f" > l'approximation de f f(t)dt fournie par Q n (f) est par défaut. 

c. La méthode est convergente, en effet si M est un majorant de \f"\ sur [a, b] on a: 

h 2 1 

|-Rn(/)| S ~r(k — a )M proportionnelle à — ~ — > 
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y=f(x) 




aY li \ 



Figure 5.5: Méthode des rectangles avec point milieu 

5.2.4 Estimation de l'erreur. 

• si /" < M sur [a, 6] on a une estimation à priori de l'erreur |j(6 — a)M 

• si /" n'est pas majorée on peut minimiser le nombre d'évaluations de /. On se limite souvent 
au cas où n est une puissance de 3 ( n: 1.3,9,27.. . . ). \Q n — Q$n\ fournit en général une 
estimation de l'erreur. 



5.2.5 Exemple 

Log2 = /j 2 f » 0.6931471806 

1 



/(*) 



/'(*) 



n*) 



\f"{x)\ <2 = M, Vu > 1 



lo 



Si n = 3 on a: h = 2 -y^ = ^, 



o = 1, 01 = 1 + l| = f , 02 = 1 +2i = | et 03 = 1 + 3| = 2 donc: 
2422. = "d'où 



1,6 6 



Qs = h E /(T«) = &(/(•*>) + /(7i) + /(•»)) = 3(7 + 9 + 11) 



478 
693 



: 0.68975468 



■ Si n = 9 on a: h = 0, Oq = 1, ai = 1 + lg 



Q1+Q2 _ 21 _ a.+a.+i 

2 18' « 2 



«2 = 1 + 2| 
i-^d'où 



1 + in et ag = 2 donc: 



?9 = fc'Ë /(-*) = M/(7o) + /(71) + ■ ■ ■ + /(-»)) = k(J| + 5? + • • • + ^^ + ■ ■ ■ + lN 



«=o 
14145553414 

20419054425 



19 21 



19 + 2i 



35 



' 0.69276237 



n 


Qn 


(6-a)^ 


\Qn-Q 3 »\ 


erreur exacte 


1 
3 
9 


0.66666667 
0.68975468 
0.69276237 


8.3 10~ 2 
4.6 ÎO" 3 
5.2 10~ 3 


2.3 10~ 2 
3.0 10~ 3 


-2.7 HT 2 
-3.4 10~ 3 
-3.9 10" 4 
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5.3 La méthode des trapèzes 

5.3.1 La formule élémentaire: 

Cette méthode consiste à remplacer la courbe de la fonction / entre a et (3 par une droite passant 

par les deux points (a,f(a)) et (/?,/(/?). Ceci est schématisé dans la figure suivante: 




Figure 5.6: Méthode du trapèze 



[^ f(t)dt = £(/(«,) + /(/?)) + R a ,g(f) 



Pour exprimer le reste R a $(f), on suppose que / a une dérivée seconde /" continue sur [a,/3] et 
introduisons la fonction auxiliaire <p(x) définie par: 



<p{x) = F(*y + x) — F(j — x) — cc(/(7 + x) + f(j — x)) — kx 7 7 : 



a + p 



(5.5) 



k est une constante choisie telle que ¥?(§) = 0. En exigeant cela on obtient pour k l'expression 
suivante: 

, _ F(/3)-F(«)-|(/(/3) + /(q)) 

~ (V2) :i 

t^(0) = (p(h/2) = 0, d'après le théorème de Rolle il existe £ e]0, /i/2[ telle que (p'fë) = 0, or 

y/(as) = -#(/' (7 + .t) - /'( 7 - x}) - Zkx 2 

Le théorème des accroissements finis appliqué à /' sur l'intervalle J7 + £,7 — £[ implique qu'il 
existe r\ ÉEJ7 + £, 7 — £[ telle que 

/'(7 + 0-/'(7-î) = 2Ç/"(»7) 

comme p'(Q = on obtient fc = -§/"(*j) et comme (p(|) = = F(/ï) - F(a) - |(/(/3) + /(a)) + 
|/"(»/)(|) 3 ) on en déduit: 

£ f(t)dt = |(/(/3) + /(a)) - g/»(„)i a < „ < /3 
c'est la fomule élémentaire pour la méthode des trapèzes. 
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5.3.2 Formule composée 

On suppose que /: [a, b\ vers IR a une dérivée seconde continue /". On considère la subdivision 

régulière de pas h de [a, b]: 

«0 = a < ai < û2 < ■■■ < a n — °i sur chaque [ai, ctj+i] on a par application de la formule 

élémentaire: 

£" +1 f(t)dt = |(/( 0j ) + /(oi+i)) - §/"fâ)i «•«<&< oj+i 
sur l'intervalle [a, o] on a donc: 

& , î=n-l ,3 i=n-l 

/ /w* = % E (/k) + /k + i)-^ E /"fê) 

= K^+^IM + ^f) -g(6-a)AA) 



r»(/3 



«»(/) 



T n (f) est la valeur approchée de f f(t)dt. R n {f) est l'erreur de méthode lorsqu'on approche 
£/(*)* par T„(/). 

La figure suivante permet de schématiser cette méthode d'approximation: 




a n =a a-| a 



Figure 5.7: Méthode des trapèzes 

5.3.3 Propriétés 

• si f = alors R n (f) = 

• si /" > l'approximation fournit par T n (f) est par excès. 

• si /" > et garde le même signe sur un intervalle [a, b] on peut encadrer l'intégrale f f(t)dt 
par: 

a„(/)< / f(t)dt<T n (f) 

• 'a 

• la méthode est convergente: la convergence est en — 



5.3.4 Estimation de l'erreur 

• si f" est Majorée par M on a: \R n (J)\ < (6 - a )^ 

• pour limiter le nombre d'évaluations de /, on peut se limiter au cas où n est une puissance 
de 2, alors dans ce cas \T^n — T n \ fournit une estimation de l'erreur ainsi qu'un test d'arrêt 
des itérrations. 
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4- Intégration et dérivation numérique 



5.3.5 Exemple 

Log2 = ff f Ri 0.6931471806 



n 


T 


ib-a)if- 


\T 2 n~T n \ 


erreur exacte 


1 


0.75 


8.3HT 2 


— 


5.710~ 2 


2 


0.7083333 


2.110 -2 


4.710" 2 


1.510~ 2 


4 


0.6970238 


5.2HT 3 


1.110" 2 


3.410~ 3 


8 


0.69412183 


1.310" 3 


2.910- 3 


9.810- 4 



5.4 Accélération de la convergence : Méthode de Simpson (1710- 
1761). 

5.4.1 Formule élémentaire 

Pour établir la formule élémentaire des rectangles et des trapèzes, on a en fait interpolé la fonction 
/: 

• par une droite horizentale ( un polynôme de degré 0) dans la méthode des rectangles avec 
point milieu. 

• par une droite ( un polynôme de degré 1 ) dans la méthode des trapèzes. 

On peut espérer améliorer la précision en approchant le graphe de / par l'arc de parabole d'axe 
vertical, passant par les points A(a, f(a)), B(f3, /(/?)) et (7(7, /(j)) (avec 7 = 2! )■ ^ e ^^ rev i en t 
à interpoler / en ces points par un plynôme de degré 2. 
On montre facilement que le poynôme de degré 2 qui passe par A : B et C est de la forme: 

P(T) (X—y)(X-l3) (X-a)(X-i3) (X - a){X - -y) 

(a-7)(a-/3) (7 - a) (7 - 0) (fi - a)(p - 7) 



l, 2 



[(X _ 7 )(x - /3)f(a) - 2(X - a)(X - 0)f{ 7 ) + (X - a)(X - 7 )/(/3)] 



On a donc f^ P(X)dX = g (/(a) + 4/(7) + /(/?)) on propose donc la formule élémentaire suivante: 



f f(t)dt = ^{f{a)+if(i)+f(fl) 



+ R a 



pour exprimer le reste, en est amené à supposer que / a une dérivé quatrième f^ A > continue: 
considérons la fonction auxiliaire: cp : [0, |] — » IR définie par 

<p{x) = F( 7 + x) - F( 7 -x)- |(/( 7 + x) + /( 7 - x) + 4/(7)) + kx b 

k est une constante choisie telle que: ip(^) = 0. 

<p{0) = fyk) — 0, d'après le théorème de Rolle, il existe Ci G [0, |] telle que <p'(t;i) = 0. 

D'après l'expression de (f'(x) donnée ci dessous 

•P'(x) = |(/( 7 + x) + /( 7 - x)) - J/( 7 ) - |(/'(x + 7 ) " fil ~ x)) + 5kx i 

on conclut que <^'(0) — 0, et comme ^'(Çi) = on en déduit du théorème de Rolle qu'il existe 
& e [0,6] teille que ip"{(, 2 ) = 



f {2 \x) = -/(7 + x) - /( 7 -x) - -(f (2) (x + 7) + / (2) (7 " x)) + 20kx 3 
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f(p) 

f(Y) 


y=f(x) 






^^^" 


/ 


f(a) 






X 




a y 


P 





Figure 5.8: Méthode de Simpson 

V^'fe) = 0, i^' 2 '(0) = 0; par application du théorème de Rolle de nouveau, il existe Ç3 € [0,^2] 
teille que <p (3) (6) = 0. 

Le même procédé appliqué à <p{> donne: 

V W( X ) = --(/< 3 >(:r. + 7) - /< 3 >( 7 - x)) + 60fcr 2 

cpP'fe) = implique |(/< 3 >(à + l) " / (3) (7 " ù)) = 60fc^ 3 2 

Le théorème des accroissements assure l'existence de r\ G] 7 — £3, 7 + ^[ tel que /' ^(£3 + 7) — 

/ m (7-&)=2fi>/ W fo) 
On a donc finalement: 

ib = JL/c*)(«) 

90 ^ w; 
avec cette valeur de A; mise dans l'exprerssion de la fonction auxiliaire (p et en demandant que 
^(2) - on en déduit: 

[ /(*)* = ^(/3) " ^(«) = J(/(0 + /(«) + 4/(7)) - ^/ (4) (>7) 
ceci est la formule élémentaire pour la méthode de Simpson. 

5.4.2 Formule composée 

On suppose que / : [a, b] vers IR a une dérivée quatrième / \ J continue. On considère une subdivision 
régulière de pas /i = ^^ de [a, 6]: ao = a < a\ < 02 < ■■■ < a n — b 
sur chaque intervalle [a^Oi+i] on a: 

£" +I /W* = *■(# - F(a) = £(/(*) + f(a l+1 ) + 4/(^+1)) - ^/ <4, fe) 
d'où sur [a, 6]: 

fj{t)dt ='f±\f(a,) + /(a 1+1 ) + 4/(^±l)) - ^' = f '/Wfe) 
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En appliquant le théorème des valeurs intermédiaire au dernier terme du membre de droite: on 

trouve: 

f\f(t)dt = ^f(a) + f(b) + 2~J2f(a,) + 4~J2f( a ^^)) - ^(b - a) *£ ' /« (0) 

i=l i=0 ' î=0 



Sn(f) 



Rn(f) 



Le reste R n (f) représente l'erreur de méthode lorsqu'on remplace / f{t)dt par la somme finie 

s„(f). 

5.4.3 Propriétés 

• Rn(î) = S h a f(t)dt-S n U) 

• la méthode est convergente, sa convergence est réglée par le terme en h 

. £„(/) = i(2Q„(/) + r„(/)) 

5.4.4 Estimation de l'erreur. 

• En pratique il est difficile d'obtenir une majoration de la dérivée 4 eme . Sinon, si M est un 
majorant de f^- 4 ' on a: 

• Pour limiter le nombre d'évaluations de /, on peut se limiter au cas où n est une puissance 
de 2. Alors dans ce cas \S-za. — S n \ fournit une estimation de l'erreur ainsi qu'un test d'arrêt 
des itérations. 

5.4.5 Exemple 

fjf = ln2 « 0.6931471806, f(x) = -^, f"(x) = J? f^(x) = $ fW(x) = f|. On remarque 
donc que f^ >(x) est majorée par 24. 



n 


nombre de points 


s„ 


(b-a)^f 


\S-2" S„\ 


erreur exacte 


1 
2 
4 
8 


3 

5 
9 

17 


0.69444443 
0.69325393 
0.69315450 
0.69314760 


8.3 10~ 3 

5.2 10~ 4 

3.3 HT 5 
2.0 10~ 6 


1.2 10~ :! 
1.1 10~" 
6.9 10~ 6 


1.2 10 -3 

1.1 10~ 4 

7.3 10~ 6 

4.2 10~ 7 



5.5 La méthode de Gauss ( à deux points ) 

5.5.1 Principe de la méthode de Gauss à deux points 

En établissant la méthode des trapèzes élémentaires: 

r h ., , . 2 „, , h 



l /(«)* = |/(«) + fc|/(7) + J/09) 



on s'est fixé les abscisses a, j3 et 7 = ^^ et on a calculé les poids §, -xh et 4 de façon à ce que la 

formule soit exacte pour les polynômes de degré inférieur ou égal à 2. 

Alors on peut se poser la question suivante: pourquoi ne pas considérer les abscisses elles-même 
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comme des inconnues à déterminer en même temps que les coefficients? 

on se limite à 2 abscisses et on propose pour approximation la formule élémentaire suivante: 

/ fl /(t)dt = A'/( 7 ') + A"/( 7 ") 

■la 

où A', A", 7' et 7" sont à déterminer pour que le reste soit nul si / est un polynôme de degré 0,1, 
2 ou 3. 

L'ensemble des polynômes de degré inférieur ou égal à trois est un espace vectoriel sur IR. Il admet 
comme base la famille B = (l,t, t , t ). Si on exige que la formule élémentaire précédente soit 
vérifiée pour les polynômes de degré inférieur ou égal à trois elle doit être vérifiée pour chaque 
élément de la base B. En éxigant cela on obtient les relations suivantes: 

(j - a = A' + A" pour f(t) = 1 
i(/3 2 - a 2 ) = A' 7 ' + A"7" Pour f(t) = t 
i(/5 3 -a 3 ) = A'7 ,2 + A"7" 2 pour f(t) = t 2 
i(/3 4 - Q 4 ) = A'7 ,3 + A"7'' 3 Pour f(t) = t 3 

pour simplifier les calculs on va établir le résultat pour a = —1, $ = 1 et avec un changement de 
variables approprié on remonte à a et /3 quelconques. On a donc le système suivant: 

2 = A' + A" (1) 

= A'V + A" 7 " (2) 

| = A' 7 ' 2 + A" 7 " 2 (3) 

_ = A' 7 ' 3 + A" 7 " 3 (4) 

La résolution de ce système permet d'obtenir la relation suivante: 

1 , „, 1 , 



| i /(t)dt = /(--^) + /( 



\/3 ' V3 
en effectuant le changement de variables u = '-H^ — "~^t : la formule précédente devient: 

on peut alors prouver la formule élémentaire de Gauss à deux points: 

pour une subdivision régulière de pas h = -^- de [a, b] si 7^ est le milieu de [a^, Qi+i]: on a donc 
la formule composée: 

une estimation de l'erreur est fournie par |C?2" — G n \. 
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5.5.2 Exemple 



flf = ln2 x 0.6931471806, f(x) = -£,, f"(x) = Jj /< 3 >(:r) = £ f^(x) = ff. On remarque 
donc que /' ' [x) est majorée par 24. 



n 


nombre de points 


G n 


(b-a)^f 


Lj2« — G n \ 


erreur exacte 


1 
2 
4 

8 


2 
4 
8 
16 


0.69230769 
0.69307663 
0.69314227 
0.69314681 


5.6 10"' ! 
3.5 10" 4 

2.2 10~ 5 

1.3 10" 6 


7.7 10" 4 
6.6 10" r ' 
4.6 10" 6 


-8.4 10" 4 

-7.1 10" 5 

-4.9 10" 6 

-4 10" 7 



5.6 Formule de quadrature 

5.6.1 Préliminaire et exemples 

Comme nous l'avons signaler au début, la plupart des algorithmes numériques procèdent comme 
suit: on subdivise [a, b] en plusieurs sous-intervalles (a = xq j x\ j x^ j - . . j x n = b) et on utilise le 
fait que: 

/ f(x)dx = J2 r +1 f^)das 

Ja j=Q Jxj 

De cette manière, on est ramené au calcul de plusieurs intégrales pour lesquelles la longueur de 
l'intevalle est relativement petite. Si on s'intéresse par exemple à: Ij = j.^ :l + 1 f(x)dx, on posons 
x = Xj + thj avec hj = Xj + i — Xj. l'intégrale Ij devient: 

Ij = / f(x)dx = hj / f(xj + jhj)dt 

Jx.j Jq 

Notons enfin que g(t) = f(Xj + jhj). Il reste alors à calculer une approximation de: 

/ g(t)dt 

■la 
Exemples: 

• Dans le cas de la formule du rectangle avec point milieu, on a 

j 9(t)dt « g(i) 

• Dans le cas de la formule du trapèze: 

/ fl(i)««4[s(0)+<7(l)] 
JO 2 

• Dans le cas de la formule de Simpson: 

j\(t)dt » ifo(0) + ig(\) + 5(1)] 
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5.6.2 Définitions 

Définition de la formule de quadrature: Une formide de quadrature à s étages est donnée 
par: 



/ g(t)dt^Y^b ig ( Ci ). (5.6) 

■h) -, 



Les Ci sont les nœnds de la formule de quadrature et les bi sont les poids. 

Exemples: 

• la formule du rectangle avec point milieu f g{t)dt ss g{\) est à un seul étage s = 1 le nœud 
Ci = h et le poids &i = 1. 

• la formule du trapèze f g(t)dt ~ 2#(^) + 2#(1) es ^ ^ deux étages s = 2. Les nœuds c\ = et 
C2 = 1, les poids fri = 62 = %. 

• la formule de Simpson: / g(t)dt sa i</(0) + is'(è) + è^CO es ^ ^ trois étages s = 3. Les nœuds 
ci = 0, C2 = 5 et C3 = 1, les poids &i = 63 = ^ et 62 = i. 

5.6.3 Ordre de la méthode 

Définition de l'ordre de la formule de quadrature: On dit que p est l'ordre de la formule 
de quadrature (5.6) s'elle est exacte pour tous les polynômes de degré p — 1, c'est à dire: 



„l s 

/ g(t)dt = y^big(Ci) pour deg# < p - 
J ° z=l 



(5.7) 



On voit que la formule du rectangle avec point milieu et du trapèze sont d'ordre 2. La formule de 
Simpson est d'ordre 4. 

Théorème: La formule de quadrature (5.6) a un ordre p si et seulement si les p conditions suivantes 
sont vérifiées: 

î> c î" 1 = - P° ur ? = 1.2,...,î> (5.8) 



1=1 



'1 



Preuve: 

— * la formule (5.7) appliquée successivement à (g(t) = t q )(q=i,2 p) donne: 



jy^-bP* 1 



La réciproque découle du fait qu'un polynôme de degré p — 1 est une combinaison de linéaire de 
1, t, t 2 ,. . . , £ p_1 c.à.d g(t) = «i + a 2 t + «:}£ 2 + • • ■ + Cï p £ p_1 et que l'intégrale f Q g(t)dt ainsi que 
l'expression Ya=1 °i9{ c i) son t linéaire en g. Par application de la formule (5.8) pour q = 1, 2, . . . ,p 

on obtient p équations. En multipliant l'équation q = 1 par CKi, l'équation q = 2 par Où2, , 

l'équation q = p par a p et en faisant la somme terme par terme on trouve: 

Qil + q 2 - + . . . + a v - = Qi / ldt + a 2 tdt + . . . + a p t v ~ x dt = 
2 p Jo ./o ^o 

= / 9(t)dt = 
■h 
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i=\ i=X 

s 



1=1 



= Y J b l g(c i ) 

i=l 

ce qu'il fallait démontrer. 

Remarque: 

Détermination des poids b{ de la formule de quadrature. 

En fixant les nœuds ci, ..., c s ÉE [0,1] (distincts), la condition (5.8) avec p = s représente un 

système linéaire pour b\ 1 . . . , b s : 



/ 1 1 

Cl c 2 

S — l S — l 

V C l C 2 



b 2 



! 



\b s ) 



(5.9) 



\\) 



La matrice dans (5.9) est une matrice de Vandermonde, les q sont distincts alors son déterminant 
est non nul. Le système admet donc une solution unique. Ceci nous donne une formule de quadra- 
ture d'ordre p = s. 



cas ou les a sont équi distant s: 

• s = 2 donc ci = et 01 = 1 on a alors le système 2x2 suivant: 



1 1 
1 



ce qui donne b\ = 62 = ô ' ^ a f° rmu l e de quadrature correspondante est la formule du trapèze: 
tig(t)dt = \g{0) + \ g (l). 

s = 3 on a C\ = 0, 02 = \ et Cg = 1, alors le système 3x3 suivant: 



1 1 
1 



iW l 2 , v -, / . , 

ce qui donne b\ = 63 = g et 62 — 5, la formule de quadrature correspondante est la formule 
de Simpson: tf g(t)dt = | s (0) + |ff(|) + Jff(l). 

s = 4 on a Ci = 0, C2 = i 5 C3 = f et C4 = 1, alors le système 4x4 suivant: 



/ 1 1 



1 \ 



ê - 1 
(i) 2 (§> (1) : 

\o (|) 3 (I) 3 (l) 3 / 



b 2 
V64/ 



/ 1\ 



v!/ 
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ce qui donne b\ = 64 = i et 62 — ^3 — §? l a formule de quadrature correspondante est la 
formule de Newton: 

jf 1 s(*)* = \g(o) + §*(J) + §$(§) + gff(i) 



On vérifie facilement que la formule de Newton est d'ordre 4. 



1, le système 5x5 correspondant 



90 ;î 90 ■ 

correspondante est la formule de Boole et est d'ordre 6. 

f 1 , x , 7 ,„, 32 ,1, 12 ,2, 32 .3, 7 ,,, 

/ g(t)dt = —o(0) + — g(-) + — g(-) + —g(-) H g(l) 

J y ^ 90 ^ ; ^9(tV 90 V 90 y V 90' V ; 



5.7 Dérivation numérique 

Dans cette partie nous abordons le problème d'approximation de la dérivée d'ordre k de / en un 
point xq donné: /' >(xq), k = 1.2.... Nous allons voir que si l'on connais les valeurs prises par / en 
un certain nombre de points Xi = Xq + ih, % = 1, 2, . . . il est possible de trouver une approximation 
de /( >(xq). Dans ce qui suit, nous supposons que / est suffisament dérivable au voisinage de xq. 

5.7.1 Approximation polynômiale 

Dans le chapitre précèdent, nous avons montré que toute fonction / : [a, b] — ► IR, / ÉE C' ïl+1 '[a, b], 
si xq G [a, b] (un point fixe), alors pour tout x G [a, 6], il existe un élément c entre xq et z tel que 

f(x) = P n (x) + E n (x) (5.10) 

iU*) = /(*o) + -^r^ - *o) + ... + ^jp^G* - *o) n (5.11) 

^^ = („ + !)! ( g - g °) ( 5 ' 12 ) 

P n est le polynôme qui peut être utilisé comme approximation de /, E n {x) est l'erreur commise 
lors de cette approximation. Si on se limite au cas ou n = 1 on aura: 

f(x) = /(.t„) + f(x )(x - x ) + ^yp (^ " *o) 2 (5.13) 

La quantité, f (xq) -\- f (xq)(x — xq) n'est rien d'autre que l'équation de la tangente de / au voisinage 
de xq. Dans ce voisinage, la droite /(xq) + J'{xq){x — xq) est une bonne approximation de f(x). 
Cette approximation est d'autant meilleure que x est très voisin de xq, c'est à dire que: Ax = x — xq 
est petit et donc l'erreur est d'ordre Ax 2 . 

5.7.2 Approximation de la dérivée première 

A partir du développement de Taylor (5.13), on est en mesure de donner une approximation de la 
dérivée première de / en ir : 

gfr)- /frD + ^-/(SD) _gg (e) (5,4) 
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L'approximation de la dérivée première est donc donnée par: 

J, s df f(x + Ax) - f(x ) 
/ + (*û) = ^(xo) « ^ (5.15) 

cette approximation est appelée différence finie progressive (où avancée). Si on remplace Ax 
par — Ax dans (5.14) on obtient: 

df ( v /(ffp - Ax) - f(x Q ) Axd 2 f 

dx M = =ÂÏ + -Tdx^ {c) (5 ' 16) 

de l'équation (5.16), on obtient l'approximation suivante de la dérivée: 

/>o) = g(so)* /( *°-** ) 3 : /( * o) (5.17) 

cette approximation est appelé différence finie rétrograde (où arrière). 

Dans les deux cas d'approximation (différence finie progressive ou rétrograde), on remarque que 

l'erreur est d'ordre 0(Ax) et est approximativement la même dans les deux cas. 

Pour obtenir une approximation plus précise de la dérivée ^:{xq), on peut faire la moyenne des 
approximations progressive et rétrograde. En effet, en faisant la somme des équations (5.14) et 
(5.16) on a: 

df f(x -Ax)-f(x -Ax) Ax t d 2 f t _^ d 2 f 
2 ^ Xo) = Âx + ^^ (C) -^ (C)1 (5 ' 18) 

puisque xq < c < xq + Ax, xq — Ax < c < xq et Ax est petit, donc c et proche de c et par 
suite la différence inh{c) — ^2 (c) tend vers zéro. On conclut que l'erreur dans (5.18) est beaucoup 
plus petite que O(Ax). Pour calculer l'erreur dans ce cas on procède comme suit: On utilise le 
développement de Taylor à l'ordre 3 de /(xq + Arc) et f(xo — Ax): 

/(x„ + Ax) = f(x„) + Ax/'(x ) + { A^.f {2) (xo) + ^^/ (3) M + ■■■ (5.19) 

/(x„ " Ax) = /(*o) - Ax/'(x ) + &*£fW( Xo ) _ i^f(3) {xo) + . . . (5 . 20) 

En faisant la différence de (5.19) et (5.20) on obtient: 

,f(x + Ax) - /(x - Ax) = 2Ax/'(x ) + ^p-/ (3) (xo) + ■ ■ • (5.21) 

On pourra montrer que (démonstration à titre d'exercice): 

/'(*„) = ti^±È±^M _ M/» (fl (5 . 22) 

ceci conduit à l'approximation suivante: 

,/, x ,/, , _ f{xo + Ax)-f{xp-Ax) 
f±(x ) = f (xq) « ^ (5.23) 

appelée différence finie centrée. Comme nous allons le voir dans un exemple explicite, la 
différence finie centrée est plus précise que la différence finie progressive et rétrograde. Ceci est 
prévisible du fait que l'erreur est d'ordre O(Ax) 2 dans le cas de la différence centrée. Dans les 
trois cas que nous avons considéré (progressive, rétrograde et centrée) le nombre d'évaluation de la 
fonction / est égal à deux. 
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Exemple 

Considérons la fonction f(x) = log(rc), on cherche à approcher f'(xo = 1) en utilisant Ax = 0.1. 
Le calcul exact de /' (1) donne /'(l) = 1. 

• A partir de l'approximation différence progressive on a: 

f+{xo = 1) fs (j{x{) + A:c = 1.1) — f(xQ = l))/(Aa:) sa 0.953102, l'erreur de cette approximation 

est de l'ordre: 

E 2 ~ |^,f (2) (l)l -0.05. 

• Dans le cas de différence arrière on a: 

f_(x Q = 1) ~ (f(x - Ax = 1.1) - f(x = l))/(-Ax) « 1.053605. 

• L'approximation de différence centrée donne: 

f ± (x = 1) « (/(aîo + &x = 1-1) - /(#0 - Aa; = 0.9))/(2A.t) « 1.00335, Terreur dans ce cas est de 

l'ordre: 

Ez ~ |^/ (3) (1)1 « 0.00333. 

Dans cet exemple, on voit clairement que la différence centrée donne une meilleure approximation 

avec une erreur de l'ordre de 10 _1 *. 

5.7.3 Approximation de la dérivée seconde 

La somme des équations (5.19) et (5.20) donne: 

f(x + Ax) + f(x - Ax) = 2/(a:o) + (Ax) 2 f"(x ) + 2 i ^-fW(x n ) + ... (5.24) 

On peut démontrer le résultat suivant (démonstration a titre d'exercice): 

fW - /('■^)-2/fa,) + /(x n -A J :) _ ÇA£ /W(G (525) 

Ceci nous conduit à la formule d'approximation de la dérivée seconde: 

/±(^o) = / (xo) ~ 7^32 ( 5 - 26 ) 

appelée différence finie centrée, cette appelation est dû au fait que l'approximation (5.26) peut être 

obtenu en appliquant deux fois la différence centrée (5.23) de la dérivée première. 

La différence centrée de la dérivée seconde nécessite trois fois l'évaluation de la fonction /: f(xo + 

Ax), f(x ) et f(xo - Ax). Les poids respectifs: , A * w , - (A ^w et ^Lp sont illustrés sur la figure 

4.9. En fait, dans toutes les formules d'approximation considérées la somme des poids doit être 

zéro. 

x -Ax x x + Ax 

• • • X 



(Ax)" 2 -2(Ax)~ 2 (Ax)" 2 Poids 

Figure 5.9: Approximation de la dérivée seconde. 



5.7.4 Approximation des dérivées partielles 

Considérons une fonction à u(x,y) deux variables x et y €E H. Supposant que / admet des dérivées 
partielles de tout ordre au voisinge de (xo,yo)- 75^ (resp ^p) n'est rien d'autre que la dérivée 
ordinaire de u par rapport à x (resp y) en gardant y (resp x constant). On pourra donc utiliser la 
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formule de différence avancée, rétrograde ou centrée pour trouver une approximation de i^ et/ou 

du 



Ox 
Ou 



(xo,yo) ' 
(xo,yo) ' 



u(x + Ax,y ) - u(x - Ax, y ) 

2Ax 
u(xg, ya + Ay) - u(x , y - Ay) 

2Ay 



(5.27) 
(5.28) 



Les deux points d'évaluation de la fonction u sont représentés sur la figure suivante: 

En physique on a souvent besoin d'évaluer le Laplacien V 2 u. A deux dimension Vu est donné 
par: 

2 d 2 u d 2 u 
dx 2 dy 2 

En utilisant la différence centrée pour 4-^ et 4~^ (5.26), on montre que V 2 u est approché par: 

u(x + Ax,y ) - 2u(x ,y ) + u(x - Ax,y ) 



V u(x ,y ) 



(Ax) 2 

u{xo,yo + Ay) - 2u(x ,y ) + u(x ,y - Ay) 

(Ay) 2 



(5.29) 



L'approximation de V 2 u nécessite cinq évaluation de u (voir figure), l'erreur de cette approximation 
est de l'ordre max(0(Ax) 2 , O(Ay) 2 ). En pratique, on utilise souvent Ax = Ay, dans ce cas V 2 u 
devient: 

V 2 u(x ,y )^ (5.30) 

u(x + Ax,y ) + u(x - Ax,y ) + u(x ,y + Ay) - iu(x ,y ) + u(x Q ,y - Ay) 

(Ax) 2 

Comme il se doit, la somme des poids relatifs est égale à zéro. 







y 








1 . 


» (0,1) 






(-1,0) 




,(0,0) 


(1,0) 


X 


1 


-4 

1< 
Poids 


» (0,-1) 


1 


P 



Figure 5.10: Approximation du Laplacien. 
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5.8 Exercices 

Exercice 1. Reprendre la formule du trapèze établie dans le cours, considérer la subdivision rég- 
ulière (ai)u—Q 2n\ ou a i — a + ^~^r- Calculer T<in en fonction de T n . Conclure que lorsque T n est 
connu, les T2 7l et T n +i nécessitent le même nombre d'évaluations de la fonction /. Montrer que le 
lien entre la formule de Simpson et celle des trapèzes est donné par 
S n = 3(4T 2n — T. n ). 

Exercice 2. On veut calculer n = 4/ rrr-s avec une incertitude absolue de 10~ 2 . On approche 
l'intégrale par les méthodes des rectangles avec point milieu, des trapèzes, de Simpson, de Gauss 
à deux points. En combien de sous-intervalles faut-il subdiviser [0, 1] dans chaque cas pour avoir 
l'incertitude désirée sur 7r? 
Faire effectivement le calcul par une méthode. On donne les renseignements suivants: 

/M = YT^ ' f{x) = - 2xf2 W 

f"(x) = 2(3.t 2 - l)f(x) , /< 3 > = 24.(1 - x 2 )f 4 (x) 

fW = 24(5z 4 - 10a: 2 + l)f(x) , / (5) = 240i:(-3:r 4 + 10k 2 - 3)f(x) 

Exercice 3 Calculer une primitive de f(t) = y/l + t 2 . En déduire I = f f(t)dt. 

Calculer / à 10~ 2 près par la méthode des rectangles avec points milieu et par la méthode des 

trapèzes. 

Améliorer la précision sur I à partir des deux approximations précédentes, on donne: 

f'(x) = x/f(x) , f"(x) = l/f(x) 

;(3) = -3x/f'(x) , /< 4 > = 3(ix 2 - l)/ f{x) 

/< 5 > = 15x(3 - 4x 2 )/f(x) 

Exercice 4. Quelle méthode proposerier-vous pour calculer: / = f () ^|x — 0.693471806| îdx à 10 -6 

près 

Exercice 5. La formule des trapèzes élémentaire pour l'approximation de l'intégrale f^ 1 f(x)dx 

par: 



r /(*)*b«J(/(*o)+/(*i)) 

Jxn 4 



Introduire la fonction auxiliaire R(h) = f^ 0+ f(x)dx — tK/(#o) + f(x$ + h)). On remarque que 



Montrer que l'erreur R = f^ 1 f(x)dx — ^{f{xo) + /(cci)) lorsque / est deux fois continuement 
dérivable a pour expression R = j^ffâ) où ^ G [xo, X\\ 
On pourra procéder de la manière suivante: 

calculer i?(0), R'(h), R"(h) et R'(0). 

On rappelle la formule de Taylor avec reste intégrale pour une fonction de lasse C 2 : q(x + t) = 

g(x) + tg'(x) + f* +t (x + t - 8)g"(0)d6. 

Exercice 6. Soit / une fonction de classe C°° sur I C IR. On considère I(f) = j f(t)dt que l'on 

désire approcher par 

J(.f) = h(af(0) + bf(^) 

1. Déterminer a et b pour que la formule approchée </(/) soit de degré de précision le plus élevé 
possible. 

2. Calculer f, 1 ^jj (h=l) 

Si f(x) = tj^m comparer /(/) et J(/). 
Exercice 7. Soit g une fonction de classe G . 
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a. Déterminer le polynôme d'interpolation de Lagrange P 2 passant par (0, #(0)), (1, <?(1)) et (2, g{2)) 

b. Trouver les constantes oùq, c*i et a 2 tels que la formule d'intégration: 



g(t)dt = a g(0) + a ig (l) + a 2 g(2) 



(5.31) 



soit exacte pour les polynômes de degré < 2. 

c. Vérifier que (1) est exacte pour P 2 

d. On posant g(t) = f(xo + jt) et en utilisant un changement de variables judicieusement choisit 
retrouver à partir de (5.31) la formule d'intégration élémentaire de Simpson sur l'intervalle [#0j#i]- 
Exercice 8. Soit la formule d'intégration suivante: 

/ f{x)dx = w f(Q) + w b f{b) , be [0,1] 

a) Déterminer wq et Wb en fonction de b pour que cette formule soit exacte pour les polynômes de 
degré < 1. 

b) Déterminer b pour que cette formule soit exacte pour les polynômes de degré < 2. 

c) Soit h € [0, 1], En déduire à partir de a) et b) que la formule d'intégration suivante: 



/' 

■Ai 



/(*) 



r/(0) ■ 



?'<!*> 



est exacte pour les polynômes de degré < 2. 

Exercice 9. On veut calculer v2 = 2 f 4 cosxdx avec une incertitude absolue de 10 _1 *. En utilisant 

la méthodes des rectangles avec point milieu, des trapèzes, de Simpson et celle de Gauss, déterminer 

le pas de la subdivision régulière qu'il faut utiliser sur [0, j] pour avoir la précision désirée. 

En appliquant la méthode de Gauss, donner la valeur de y2 à 10~ 3 . 

Exercice 10: Méthodes des rectangles à gauche. 

Soit / : [a, b] — > IR telle que /' existe et soit continue; soit a < j3 ans [a, b], 

a. Quelle est la formule d'intégration élémentaire suggérée par la figure (a)? 



v=l(x) 



«B 

f(a) 



_fc_ -_ÉL. 



Figure 5.11: Méthodes des rectangles à gauche figure (a) et à droite figure (b). 

Quelle est l'erreur élémentaire? 

b. Soit ûq — a < a i < a 2 < ■ • ■ < a n — b une subdivision régulière de pas h = ^^ de [a, 6]. Quelle 
est la formule d'intégration composée obtenue en répétant la formule simple sur chaque intervalle 
[ai,aj+i]? Erreur de méthode? Que dire de l'approximation si / est monotone sur [a, 6]? Justifier 
vôtre réponse. 

c. Application au calcul de J^ ~. Calculer n pour que l'erreur de méthode soit inférieur à 10~ 3 . 

d. Reprendre les mêmes questions dans le cas de la figure (b). 



Exercice 11. 



^ETtfUP 



4. 6 Exercices 85 

On désire calculer ln2 = f Q ^r avec une incertitude absolue de 10~ 2 . 

En utilisant la méthode des rectangles avec point milieu, des trapèzes et de Simpson déterminer le 

pas de la subdivision régulière qu'il faut utiliser sur [0, 1] pour avoir la précision désirée. 

En appliquant la méthode de Simpson, donner la valeur de ln2 à 10 -2 . 

Exercice 12. 

Soit / une fonction de classe C 2 [a, b], on pose 7 = ^2~- 

1. Ecrire le polynôme d'interpolation de Lagrange P^it) passant par les points Aq = (a, /(a)), 
Ai = hJ(f))etA 2 = (b,f(b)). 

2. Montrer que: 



■'a 



-" = ^[/(«») + 4/(7) + /(»)] 



3. On considère la formule d'intégration suivante: 

j_ f(t)dt = Oi/(-l) + 03/(0) + 03/(1) + «(/) (5.32) 

R(f) étant le reste. Déterminer ai, a% et 03 pour que la fomule d'intégration (5.32) soit exacte 
pour les polynômes de degrés < 2. 

4. Pour les valeurs de ai, a 2 et a% trouvées dans 3, la formule d'intégration (5.32) est elle exacte 
pour les polynômes de degrés < 3. 

5. A l'aide d'un changement de variable, retrouver à partir de (5.32) la formule de Simpson: 

" '"[/(«) + 4/(7) +/(6)] 



/ f(t)dt: 
■ 'a 



6 

Exercice 13. 

Soit / une fonction de classe C 2 sur [a, b], soit [a, (3\ un intervalle élémentaire contenu dans [a, b]. 

a) Déterminer le polynôme de Lagrange P(X) interpolant / aux points d'abscisses a et (3. 

b) Donner et sans le démontrer l'expression de f(x) — P{x) en fonction de x, a, et /". 

c) En déduire de b) la formule élémentaire d'intégration des trapèzes. 

d) Soit (a,:)i=o,...,n une subdivision régulière de [a, b] de pas h = (b — a)/n. Retrouver la formule 
des trapèzes composée à partir de c). 

Exercice 14: 



x < 1 

^0 



F(x)= f g(t)dt, 
Jq 

a.) Quelle est la valeur de F(x) en x = | 

b.) Donner le degré et l'expression du Polynôme de Lagrange qui interpole g(x) aux points 0, 3, 

2 
3 - 

c.) Trouver des coefficients cq. c\ et C2 tels que pour tout polynôme P de degré inférieur ou égal 
à 2, on ait 



r 

J Q 



P{x)dx = c P(0) + ciP(l) + c 2 P(2) 



cl.) En utilisant un changement de variable, déduire de la formule précédente les coefficients do, 
d\ et di tels que pour tout polynôme Q de degré inférieur ou égal à 2, on ait 



/ 3 Q(x)dx = d Q(0) + d x Q{\) + d 2 Q(l) 
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86 4- Intégration numérique 

e.) Utiliser cette formule pour donner une valeur approchée de ln3. 
Quel est l'ordre de cette méthode? 



ETIWJR 



com 



eu 



eu 



Programmation D" m .. '^ 

yi p qj flJ Algèbre 2 

Coursai e $ $ " 

:*S Résumés .pS6 "" -" 



O -r-. 



D" > 



<D 



CU 



o_* 



■^.O" Analyse c Diapo u ^ Ë 

tX6iCIC6S a; g- 

Contrôles Continus ■£ ^ 

Langues mtu^ orf 

Thermodynamique — ^, ^ S 

Multimedia UlVCTS 

Economie Travaux Dirigés ±i 

Chimie Orqanîque 2 



et encore plus» 



